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(A1) Translate this text EEG leads are connected to a patient's head (14) by a set 
of surface electrodes which transmit signals over a patient cable (16) to a 
19-channel EEG data acquisition system (12). Data acquisition system (12) 
filters, amplifies and digitizes the EEG waveforms and sends the digitized data to 
the microcomputer (18) via high speed synchronous serial line (26). 
Microcomputer (18) processes the serial data stream to generate computed data 
arrays. These arrays are then used in conjunction with predetermined reference 
arrays derived from clinical studies to produce diagnostic indices. These indices 
are displayed on graphics display (20). Printed output of the diagnostic index is 
available on the hard copy output device (22) which is connected to the 
microcomputer (18). Interaction between the operator and the system is provided 
by means of a keyboard (24) and a pointing device (28) with feedback provided 
by graphics display (20). 
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CEREBRAL BIOPOTENTIAL ANA LYSIS SYSTEM AND METHOD 
CROSS REFERENCE TO RELATED APPLICATION 

This application is a continuation-in-part application 
of co-pending United States Serial No. 107,357 filed Oc- 
tober 9, 1987, which is also assigned to the assignee of 
the present application 

BACKGROUND OF THE INVENTION 

The present invention relates to a real-time, high-fre- 
quency, high-resolution cerebral biopotential analysis sys- 
tem and method, and more particularly to a computer-based 
biopotential diagnostic system and method for quantitative- 
ly determining, in a noninvasive manner, cerebral phenomena 
that can be ascertained by analyzing cerebral electrical ac- 
tivity. 

Despite a considerable expenditure of time and effort, 
current approaches to the quantitative, noninvasive assess- 
ment of cerebral electrical activity, as displayed in an 
"EEG" waveform, have not been successful in fully extract- 
ing all of the information which is present in this complex 
waveform. A great need remains for an accurate, sensitive, 
reliable, and practical neurologic profiling technology. In 
particular, contemporary intra-operative EEG monitoring 
techniques have not been widely adopted due to their in- 
herent limitations. Indeed eighty per cent (80%) of all 
medical malpractice suits are believed to be related to 
post-anesthesia morbidity and mortality, and if such EEG 
monitoring techniques were reliable they certainly would 
have been adopted. 
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A number of devices known in the prior art are capable 
of tracking cerebral activity qualitatively. Techniques in- 
volving the use of the "classical", conventional analog EEG 
are restricted to analyses in the time domain, and require 
considerable training for adequate interpretation. 
Moreover, since the frequency resolution of the human eye 
at standard speeds and gain is 30 - 60 Rz, much high fre- 
quency content is invisible. Thus visual EEG assessment is 
better characterized as being an art rather than a science. 
In fact, it has been shown that the average correlation be- 
tween seven experienced readers did not exceed 56 per cent. 

The use of frequency (power spectrum) analysis of the 
EEG in the 1960's introduced the notion of some basic 
processing of the signal prior to visual inspection and led 
to the application of frequency analysis of the EEG to 
various cerebral monitoring problems. In the past 25 years 
at least 100 papers have been published in the medical 
literature describing applications of power spectral 
analysis for purposes such as assessing the depth of anes- 
thesia and cerebral ischemia under various intraoperative 
conditions. United States Patent No. 4,557,270 issued to 
John also describes the use of power spectral analysis to 
evaluate cerebral perfusion during open heart surgery. 
Several recent studies, however, have shown many deficien- 
cies in the use of power spectral analysis to monitor 
cerebral perfusion and to determine post operative 
neurologic outcome. In addition, neither power spectral 
analysis nor any other monitoring technique has been shown 
to be reliable, and this is demonstrated by the fact that 
the well-accepted Harvard Medical School Anesthesia Monitor- 
ing Standard does not include any type of intraoperative 
neurologic monitoring due, in all likelihood, to the com- 
plexity of interpreting raw EEG data and the unreliability 
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of existing automated systems utilizing power spectral or 
time-domain analytic techniques. 

The discharge of thousands of bio-electrically active 
cells in the brain, organized in larger, interacting neural 
centers contributes to the formation of an electrical sig- 
nal with a wide frequency spectrum and extremely complex 
dynamics. Embedded in that signal is information regarding 
frequency content, non-linearities, and phase relationships 
arising from the complex neuronal firing patterns that take 
place. Because of the complexity of the EEG signal, conven- 
tional time and frequency modes of analysis have not been 
adequate to fully profile its behavior. In the Fourier 
transform of the second order autocorrelation function (the 
power spectrum) processes are represented as a linear summa- 
tion of statistically uncorrelated sine-shaped wave com- 
ponents . Contemporary approaches to monitoring the EEG by 
means of the power spectrum have thus suppressed informa- 
tion regarding non-linearities and inter-frequency phase 
relationships and are of limited utility in representing 
the EEG' s dynamic structure. Furthermore the high frequency 
low amplitude elements of the EEG have been discarded to 
date by the filtering and sampling characteristics of known 
analysis techniques. 



Because the EEG has a wide spectrum and is highly 
dynamic and non-linear, the phase relationships within the 
EEG, especially in the higher frequencies, must carry a 
great deal of diagnostic information regarding cerebral 
function. The Fourier transform of the third order autocor- 
relation function, or autobi spectrum, is an analytic 
process that quantifies deviation from normality, quadratic 
non-linearities and inter- frequency phase relationships 
within a signal. The Fourier transform of the third order 
crosscorrelation function, or crossbispectrum, is an 
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analytic process that provides similar information between 
two signals. 

Aut obi spectral analytic techniques have been applied 
to the EEG signal and the basic bispectral properties of 
the conventional EE6 focusing on frequencies below 32 Hz 
have been investigated. Such studies have also been con- 
ducted to search for changes between waking and sleeping by 
means of autobispectral analysis. Autobispectral analysis 
and power spectral analysis have also been used in an at- 
tempt to show that the EEGs of monozygotic twins are 
similar in structure. 

To date, ©o previous study has examined the high fre- 
quency (greater than 32 Hz) content of the EE6 and found in- 
formation of diagnostic value. It also does not appear that 
any study has shown autobispectral or crossbispectral 
analysis to be of any value for any diagnostic purpose and 
certainly neither of these analytic techniques have been 
shown to have any value in quantifying depth and adequacy 
of anesthesia, pain responses induced by surgical stress, 
cerebral ischemia, consciousness, degrees of intoxication, 
ongoing cognitive processes or interhemispheric dynamic 
phase relations. 

It is therefore a principal object of the present in- 
vention to provide a noninvasive high resolution high fre- 
quency electro-encephalographic system and method capable 
of recognizing and monitoring physical phenomena that are 
reflected in cerebral electrical activity. 

Another object of the present invention is to provide 
a noninvasive electroencephalograph^ system and method 
capable of determining and monitoring depth and adequacy of 
anesthesia, pain responses during surgical stress, acute 
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cerebral ischemia, level of consciousness, degrees of in- 
toxication and normal or abnormal cognitive processes. 



SUMMARY OF THE INVFNTIOM 



Accordingly, the system and method of the present in- 
vention utilises a suitable electrode and amplifier system 
to obtain 19 unipolar EEG signals from regions of interest 
on both left and right hemispheres of a subject's brain. 
Band-pass filtering of 2 - 500 Hz is used to obtain signals 
with a high frequency content. High gain amplifiers maxi- 
mize the dynamic range for the high frequency, low energy 
wave components of the signals. The system applies digital 
sampling techniques to the signals and transmits digitized 
data over a high speed serial line to a host computer. The 
system divides a 32 second long data segment from each lead 
into 128 consecutive 0.25 second intervals. The system nor- 
malizes all 19 unipolar leads by the standard deviation, 
and then characterizes the dynamic phase relations within 
the signal by processing for autobispectral variables using 
either a Fast Fourier Transform (FFT) based approach, or a 
parametric cubic fitting approach. Similarly three cor- 
responding left and right hemisphere data pairs are normal- 
ized in the same manner and dynamic phase relations between 
two hemispheres are then characterized by processing for 
crossbi spectral estimates utilizing either the FFT or 
parametric based techniques. The outcome is a set of two 
dimensional arrays representing the dynamic interactions be- 
tween all the possible combinations of frequencies (frequen- 
cy pairs) in the spectrum of interest. For each unipolar 
lead, three arrays are produced: autobicoherence, 
autobispectral density and autobiphase. Three arrays are 
also generated for each bipolar data set: crossbi coherence, 
crossbispectral density and crossbiphase . 
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Each of the autobi spectral and crossbi spectral arrays 
contains 16,512 data points. Although all, or nearly all, 
of these values can be expected to change from normal 
during different interventions or due to differing disease 
states, in the preferred embodiment only those points which 
show the greatest fidelity in tracking the particular diag- 
nostic determination in question are utilized to create a 
diagnostic criterion. The ensemble of points most sensitive 
to a particular intervention or ongoing physiologic process 
can be used to create a clinically useful single-number 
index from the computed bispectral arrays. The system uses 
these indices as a diagnostic figure of merit for the as- 
sessment of depth and adequacy of anesthesia, pain respon- 
ses during surgical stress, acute cerebral ischemia, level 
of consciousness, degree of intoxication and normal or ab- 
normal cognitive processes. This approach makes it possible 
for any, even unskilled, operator to meaningfully interpret 
the output ^of the diagnostic device. 

In situations where continuous monitoring is required, 
indices can be continuously displayed on a video terminal 
thereby enabling the operator to interactively evaluate 
regions of interest. For record keeping purposes index 
values and other pertinent variables can be sent to a hard 
copy output device or stored on a disk. 

These and other objects and features of the present in- 
vention will be more fully understood from the following 
detailed description which should be read in light of the 
accompanying drawings in which corresponding reference 
numerals refer to corresponding parts throughout the 
several views. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a schematic view of the system of the 
present invention for detecting cerebral phenomena in a non- 
invasive manner; 

Fig. 2 is a schematic view of a 19 channel EEG data ac- 
quisition system including a serial interface utilized in 
the system of Fig, 1; 

Fig. 3 is a schematic view of the microcomputer used 
to calculate and display the EEG bispectrum in the system 
of Fig. 1; 

Fig. 4 is a schematic view of the processing opera- 
tions performed by the system of Fig. 1; 

Fig. 5 is a flow chart of the operations of the 
monitor module shown in Fig. 4; 

Fig. 6 is a view of a sample display representation of 
bispectral values generated by the system of Fig. 1; 

Fig.. 7 is a flow chart of the operations of the ac- 
quisition and EEG raw data management module of the system 
shown in Fig. 4; 

Fig. 8 is a flow chart of frequency domain based 
method for producing autobi spectrum or crossbispectrum used 
by the system of Fig. 1; 

Fig. 9 is a flow chart of a parametric based method 
for producing autobi spectrum or crossbispectrum in the sys- 
tem of Fig. 1; 
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Fig. 10(a) is an illustration of a graph showing a 
bi spectral density array generated by the system of Fig. l; 

Fig. 10(b) is an illustration of a graph showing a 
biphase array generated by the system of Fig. 1; 

Fig. 10(c) is an illustration of a graph showing a 
bicoherence array generated by the system of Fig. 1; 

Fig. 11 is a flow chart of the diagnostic index genera- 
tion module shown in Fig. 4; 

Figs. 12(a) - 12(c) are illustrations of arrays of 
bispectral density values for three different states of one 
patient; 

Figs. 13(a) - 13(b) are graphs of statistical arrays 
generated by the system and method of the present invention; 

Fig. 14 is an annotated continuous aut obi spectral den- 
sity diagnostic index graph for one lead generated by the 
system of Fig. 1. 
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DETAILED DESCRIPTION OFTHF PR FFERRED EMROn iMPMTQ 

Referring to Fig. l the apparatus of the present inven 
tion includes a 19 channel EEG data acquisition system 12 
connected to a microcomputer 18 through a high speed serial 
interface 26. 



The EEG leads are connected to a patient's head 14 by 
a set of surface electrodes. The International 10/20 
electrode system and nomenclature is preferred. The EEG sig- 
nals are picked up by the electrodes and transmitted over a 
patient cable 16 to the EEG data acquisition system 12. 

The data acquisition system 12 filters, amplifies and 
digitizes the EEG waveforms and sends the digitized data to 
the microcomputer 18 via a high speed synchronous serial 
line 26. In addition, the serial line 26 can be used to 
download filtering, gain and sampling rate instructions to 
the data acquisition unit 12. 

The microcomputer 18 processes the serial data stream 
in order to generate all computed data arrays. These arrays 
are then used in conjunction with predetermined reference 
arrays derived from clinical studies to produce diagnostic 
indices which indicate the status of the patient. These in- 
dices are displayed on the graphics display 20. Printed out- 
put of the diagnostic index is also available on the hard 
copy output device 22 which is connected to the microcom- 
puter 18. Interaction between the operator and the acquisi- 
tion and analysis components of the system is provided by 
means of a keyboard 24 and pointing device 28 with feedback 
on the graphics display' 20. 

The 19 channel data acquisition system 12 is shown in 
greater detail in Pig. 2. The EEG surface potential, 
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detected by a surface electrode mounted on the patient head 
14, passes through an electrosurgery protection circuit 30, 
a defibrillator protection circuit 32, and an amplifier/ 
filter circuit 36 before being passed on to the multi-chan- 
nel analog to digital converter 38. 

The electrosurgery protection circuit 30 includes a 
radio frequency (rf) filter, which limits the rf current 
through the patient leads 16 to less than 100 microamperes 
and thus protects the patient 15 from rf burns and protects 
the amplifiers 36 from damage resulting from exceeding the 
absolute maximum input voltage specified by the manufac- 
turer. This circuit can be an LC section circuit consist- 
ing of a generic inductor connected in series to a generic 
capacitor which is then connected to ground. 

The defibrillator protection circuit 32 limits the 
voltage to the amplifiers 36 to a safe level when a 
defibrillator is applied to the patient 15 and discharged. 
This circuit can consist of a neon light bulb and or a 
parallel variable resistor connected in series to a 
grounded resistor. 

The amplifier/ filter circuitry 36 is controlled by 
the microprocessor 34 for default gain and filtering levels 
or alternate gain and filtering levels as requested by the 
operator. Preferred gain and filtering settings are dis- 
cussed later. This circuit section consists of three 
stages; the first is a pre-amplifier stage that can be as- 
sembled uiing a wide variety of high impedance pre- 
amplifiers such as those sold by National Semiconductor, 
Sunnyvale CX; the second is a programmable filters stage 
which can utilize components from Frequency Devices, Haver- 
hill MA; the third stage is a programmable amplifiers stage 
which can be assembled from operational amplifiers used in 
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conjunction with a multiplying digital to analog (D/A) con- 
verter both components can be supplied by National Semicon- 
ductor. The multiplying D/A is used to set the gain to the 
appropriate levels requested by the microprocessor 34. 

The high impedance pre-amplifier of each channel will 
saturate to either the positive or negative supply voltage 
if the input of the pre-amplifier is not terminated. This 
will lead to large positive value or a large negative value 
at the output of amplifier/ filter section 36. Such values 
will be used to identify lead failure. 

The output of all 19 channels of the amplifier/ filter 
36 is fed to a multi-channel analog to digital converter 
(A/D) 38 which is under microprocessor 34 control for sam- 
pling rate settings. The analog signals are converted to 
digital data format suitable for input to a computer. A/D 
converters sold by Analog Devices, Norwood MA can be used 
for this purpose. 

The multi-channel A/D converter 38 is optically 
coupled to data bus 40 by optical isolator 42. All control 
lines to the sample and hold circuits, the multiplexer and 
the A/D convertor 38 are also optically isolated by optical 
isolator 44. Any known optical isolators can be used for 
this purpose. 

All DC power lines going to the amplifiers 36, sample 
and hold circuits, multiplexer and A/D convertor 38 are 
also isolated from the AC power line with a DC/DC convertor 
46 in order to provide complete patient isolation from 
ground. DC/DC converters available from Burr Brown can be 
used for this purpose. 
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The basic instructions for controlling operation of 
the microprocessor 34 are stored in a read only memory 
(ROM) 48. The random access memory (RAM) 50 is used as a 
buffer memory for data and a portion of the RAM 50 can also 
be used as program memory when a control program is being 
downloaded from the microcomputer 18. 

Serial interface 52 operates under the control of the 
microprocessor 34 . The serial interface 52 is optically 
coupled with optical isolators 54 to high speed synchronous 
serial drivers 56 to provide a synchronous serial link be- 
tween the 20 channel data acquisition system 12 and any com- 
patible high speed synchronous serial interface card on any 
computer. The serial lines are isolated by optical 
isolators 54 and DC/DC convertor 58 to provide increased 
patient safety and to protect the host computer 18 from any 
transients. 

The host or microcomputer 18 of Pig. 1 is shown in 
greater detail in Fig. 3. The entire microcomputer system 
runs under control of a microprocessor 62 with the program 
memory for the microprocessor 62 being stored in ROM 64. 
The RAM 66 is used for storage of intermediate data. The 
mass storage device 84 is used for storing clinical 
databases as well as archiving patient data. 

Zn a preferred embodiment, the microcomputer 18 con- 
tains an array processor 68 (such as the Vortex sold by SKY 
of Lowell, MA) on which complex arithmetic calculations can 
be performed on entire arrays of data simultaneously. The 
preferred embodiment also includes a math coprocessor 70 
which is connected directly to microprocessor 62. The math 
coprocessor 70 is used for scalar and graphic calculations 
while the array processor 68 is used to calculate 
bispectral and other data vectors . 
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A graphics controller 72 operating under program con- 
trol of the microprocessor 62 drives a graphics display 20. 
A keyboard controller 74 interfaces directly with the 
operator's keyboard 24. A serial port 80 interfaces with a 
pointing device 82. 

Operator control of the entire acquisition, analysis 
and display procedure is controlled by the keyboard 24 and 
pointing device 82 with feedback on the graphics display 
20. One high speed synchronous serial port 76 is provided 
to interface with the 20 channel data acquisition system 
12. Port 76 can be used to send control data to the system 
(e.g., filtering, gain, sampling rate, start/ stop acquisi- 
tion, perform self diagnostics) and to receive EEG data 
from the system, as well as to download program data to the 
system. Another serial or parallel port 78 is provided to 
drive a hard copy output device 22 for printing desired 
diagnostic indices. 

Referring now to Fig. 4, a block diagram of the system 
operations and the method of the present invention is 
described. As mentioned above, the system and method of the 
present invention computes dynamic phase and density rela- 
tions of EEG signals from a preselected number of leads (19 
unipolar and 6 bipolar in the described embodiment) . Single 
number diagnostic indices are then generated from the data 
arrays by utilizing predetermined reference arrays. The 
results are quantitative indices useful for analyzing 
cerebral electrical activity as it relates to, for example, 
the assessment of depth and adequacy of anesthesia, pain 
responses during surgical stress, acute and chronic 
cerebral ischemia, level of consciousness, degree of 
cerebral intoxication and normal or abnormal cognitive 
processes . 
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Th* monitor module 402, handles the overall operations 
of the system via integration of data and process informa- 
tion from the user interface module 404, acquisition and 
raw EEG data management module 406, bispectral processing 
module 408 and diagnostic index derivation module 410. A 
detailed illustration of module 402 can be found in Pig 5. 

The user interface and display management module 404 
represents the means through which the operator controls 
and interacts with the system during the course of a proce- 
dure. This includes, but is not limited to, entry of infor- 
mation regarding the patient, type of diagnostic procedure 
being carried out, lead and acquisition settings; con- 
tinuous display of acquisition status, lead integrity, and 
diagnostic indices corresponding to regions probed by each 
electrode; and requests for printing and archiving results 
to disk. Module 404 directly interacts with the monitor 
module 402. The operations handled by module 404 can be 
achieved under one of many commercially available environ- 
ments such as Microsoft's Windows. 

The acquisition and raw EEG data management module 
406, handles all of the raw EEG data checking and process- 
ing prior to bispectral analysis. This includes, but is not 
limited to, continuous acquisition of EEG data and the 
verification of the integrity of the data; preparing all 
unipolar EEG data for aut obi spectral processing; preparing 
all bipolar EEG data for crossbi spectral processing. Module 
406 directly interacts with the monitor module 402. A more 
detailed description of module 406 is provided below in con- 
nection with Pig 7. 

The bispectral processing module 408 controls the 
generation of all data arrays measuring dynamic phase and 
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density relations within the EBG. This information can be 
organized in both autobi spectral and erossbi spectral arrays 
utilizing either an PPT based or parametric based approach. 
The tas)cs performed by this module include, but are not 
limited to: Fourier transformation; and the generation of 
power spectrum, autobispectral density, erossbi spectral den- 
sity, autobicoherence, crossbicoherence, autobiphase, and 
crossbiphase. Module 408 directly interacts with the 
monitor module 402, and a more detailed description of 
module 408 is provided below in connection with Pigs. 8 and 
9. 

The diagnostic index derivation module 410 generates 
the data values utilized in the diagnostic process. The 
task includes, but is not limited to, identifying frequency 
pairs of interest through the use of predetermined clinical 
reference arrays and creating a diagnostic index from the 
values in the bispectral data arrays at the frequency loca- 
tions defined by the reference array. Module 410 directly 
interacts with the monitor module 402, and a more detailed 
description of module 410 is provided below in connection 
with Fig. 11. 

Referring now to Fig. 5, the operation of the monitor 
module 402 will now be discussed. In step 502, the data ar- 
rays used to store the digitized EEG, the 128 0.25 second 
EEG data records, and the bispectral data of each lead are 
initialized. The data files required for storage and files 
containing data bases required for the computation of diag- 
nostic indices are also opened in the initializing step 502. 

In step 504 the system requests the information re- 
quired to start the acquisition and diagnostic process from 
the user via the user interface module 404. This requested 
information includes patient descriptive statistics (sex, 
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age, clinical symptoms ate..)/ typa of diagnostic procedure 
to be conducted, and the leads usea for aut obi spectral 
analysis and the leads used for crossbispectral analysis. 

The system includes a default mode of operation and in 
this default mode the system continuously monitors the 
depth and adequacy of anesthesia, and any pain responses 
during surgical stress utilising a default autobispectral 
density database. Default band pass filtering is performed 
from 2 to 500 Hz; the default sampling rate is set at 2000 
Hz; and default gain is automatically adjusted to achieve 
maximum dynamic range in each lead. The following discus- 
sion of the monitor module 402 will utilize the default set- 
tings of the system. 

The EEG signals measured by leads Fpl, Pp2, F7, F3, 
Fz, F4, F8, T3, C3, Cz, C4, T4, T5, P3, Pz, P4, T6, 01, and 
02 (Al or A2 f or "reference) are used for autobispectral 
analysis. 

The EEG signals measured from the difference of leads 
F7 and T3 (F7-T3) and the difference of leads F8 and T4 (F8- 
T4) originate from the area covered by the frontal left 
hemisphere 1 and frontal right hemisphere regions respective- 
ly. These signals from F7-T3 and F8-T4 are paired and used 
for crossbispectral analysis. In this way, the interhemis- 
pheric relationships for the frontal region can be ex- 
amined. Similarly, pairing C3-Cz with C4-Cz and T3-T5 with 
T4-T6 for crossbispectral analysis purposes allows for the 
examination of the interhemispheric relationships of the oc- 
cipital and parietal regions respectively. 

In step 506, 128 0.25 second buffers of artifact free 
raw EEG data are acquired. All channels transmitting ar- 
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tifactual data are properly signaled to the operator to cor- 
rect the problem. 

The system, in step 508, computes autobispectral ar- 
rays for leads Fpl, Fp2, P7, F3, Fz, F4, F8, T3, C3, Cz, 
C4, T4, T5, P3, Pz, P4, T6, 01, 02, and crossbispectral ar- 
rays for leads F7-T3 paired with F8-T4, T3-T5 paired with 
T4-T6, and C3-Cz paired with C4-Cz. Other leads may, of 
course, be used in the computation of these arrays, and two 
different approaches for bispectral computation will be dis- 
cussed below with reference to Figs. 8 and 9. 

In step 510, the single number diagnostic indices from 
all generated autobispectral and crossbispectral arrays are 
computed. Autobispectral density and crossbispectral den- 
sity clinical reference arrays are utilized in these diag- 
nostic index computations. The generation of the reference 
arrays is discussed later. The system instantaneously dis- 
plays, in step 512, all computed diagnostic indices for all 
leads being analyzed. In step 514, the system checks for an 
exit request, and if such a request has not been made, the 
system, in step 516, acquires a new 0.25 second buffer and 
repeats steps 508 through 514. in step 518, requested prin- 
touts are produced, results are stored to disk for archival 
purposes and all files are closed. In step 520, the process 
is terminated. 

A sample display representation generated by the sys- 
tem is shown in Fig. 6. Representations of the patient's 
head are shown on the graphics display in Fig 6(a) and Fig. 
6(b). The first illustration Fig. 6(a) is divided into 
nineteen sections each representing the region probed by an 
electrode. The second illustration Fig. 6(b) is divided 
into three horizontal sections representing combined left 
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and right hemisphere activity probed by the group of 
electrodes in that region. 

For head representation Fig. 6(a), each section con- 
tains a compressed continuous tracing 602 of the computed 
diagnostic index utilizing the unipolar EEG data acquired 
from the electrode in that site. For head representation 
Fig. 6(b), each section contains a compressed continuous 
tracing 604 of the computed crossbi spectral diagnostic 
index utilizing bipolar EEG data acquired from several 
electrodes in that site. 

At the request of the operator any site can be dis- 
played as an enlarged view 606 for closer examination. The 
background of the tracing of each site (such as 602 or 604) 
is color coded to reflect the possible values allowed for 
in the range of the selected diagnostic index. The most 
current value of the diagnostic index for that site will 
dictate what color is displayed in the background (e.g. Red 
= lowest value to Green « highest value) . This will 
facilitate the examination of the patient's status at a dis 
tance. Each site will be covered by a large "X" 608 if a 
lead fail or an artifact was detected for any of the leads 
contributing to the data required to generate the diagnos- 
tic index at that site. 

Referring to Fig. 7, the acquisition and raw EEG data 
management module 406 will now be described in greater 
detail. In step 702, the system checks whether the 0.25 
second buffer for which data is to be acquired is the first 
buffer being filled for that run, and if it is, the acquisi 
tion system 12 in step 704 is supplied with requested fil- 
tering, gain, sampling rate and lead selection information. 
The default settings are band pass 2 - 500 Hz for filter- 
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ing, 50,000 for gain, 2,000 samples/sec for sampling rata 
and signals from all 19 leads ara acquired. 

In step 706, the acquisition system 12 acquires data 
for each 0.25 second buffer for all requested leads and 
transfers this data to the host computer 18. The system 
detects lead fails during the acquisition cycle in step 708 
by checking for very large positive or negative values. 
Also in step 708 a publicly available algorithm is used to 
check for artifact in each lead. In step 710, leads 
generating failed and artifactual data are marked for the 
monitor module 402. 

In step 712, the system normalizes the records of data 
acquired from all artifact free leads by subtracting the 
mean of the samples in each record from each sample in that 
record, and then dividing the sample by the standard devia- 
tion of the records. This normalization sets the variance 
in each record to 1 and has the effect of weighing each 
record equally during bi spectral averaging. The process is 
therefore less dependent on the absolute power spectral den- 
sity at any frequency band. 

In step 714, each 0.25 second record from each of the 
leads Ppl, Pp2, F7, F3, Pz, P4, P8, T3, C3, Cz, C4, T4, T5, 
P3, Pz, P4, T6, 01, 02 is assigned to an Xi(t), where Xi(t) 
are the individual time series records provided for 
autobispectral processing. Also in step 714, the frontal 
left hemisphere time series, Xi(t), from F7-T3 and the fron- 
tal right hemisphere time series, Yi(t), from F8-T4 are 
provided for crossbi spectral processing. Similarly, by pair- 
ing leads C3-Cz with C4-Cz and T3-T5 with T4-T6, the cross 
bispectrum of the left and right occipital and left and 
right parietal regions respectively can be processed. It 
should be noted that for autobispectral analysis Yi(t) is 
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set to gqual Xi(t) and in all cases the index i denotes the 
record number from 1 to 128. 

In step 716, a circular buffer mechanism is used for 
storing the appropriate Xi(t) and Yi(t) records for each 
lead. The buffer is updated by storing the most recently 
acquired data record in the location of the least recently 
acquired data record. In step 718, the program checks 
whether the circular buffer has 128 acquired data records 
to start bispectral analysis, and if there are 128 data 
records in the buffer, operation of the system returns to 
the monitor module 402 in step 720. 

Referring now to Fig. 8, the frequency domain based 
procedures for producing the aut obi spectrum or the 
crossbispectrum will now be discussed. In step 802, the sys- 
tem checks whether the computation to be performed is an 
autobispectral or crossbispectral computation. 
Aut obi spectral analysis is a special case of 
crossbispectral analysis and therefore different rules of 
symmetry apply . 

In step 804, the system sets the following symmetries 
in order to proceed with autobispectral computation: 

f 1 + f2 < N/2 

where N » 512 (0.25 sees * 2000 samples in a preferred em- 
bodiment) , an d 



0 < f 2 < f l 
Xi(t) « Yi(t) > Xi(f) = Yi(f) 
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whera £1 and £2 (also referred to as Fi and F2 or Frequency 
1 and Frequency 2) denote the frequency pairs over which 
bispectral computation will be carried out. Xj.(t) and Yi(t) 
denote the individual time aeries records used for 
bispectral computation. Xi(f) and Yi(f) denote the Fourier 
transform of the time series records and i denotes the 
record number and in this embodiment ranges from 1 to 128. 

In step 806, the following symmetries are adhered to 
for crossbi spectral analysis: 

f 1 + £2 < N/2 
0 < f 1 < N/2 
0 < f 2 < N/2 
-2f 2 < £1 
Xi(t) * Yi(t) > Xi(f) * Yi(f) 

where all variables represent the same values as they do 
for autobiapectral analyaia, except that for 
crossbispectral analysis Xi(t) and Yj.(t) represent in- 
dividually derived time series records from left and right 
hemisphere leads respectively. 

The fast Fourier transform (FFT) Xi(f) and Yi(f) of 
each record of the 128 selected records for that lead, is 
computed using a standard IEEE library routine or any other 
publicly available routine in step 808. 
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In Step 810, the power spectra Pxi(f) and Pyi(f) of 
each record of the 128 selected records for that lead is 
computed by squaring the magnitudes of each element of the 
Fourier transform Xi(f) and Yi(f) respectively. 

The system computes the average complex triple product 
in step 812 by utilizing the following equations where 
bCi(fi,f2) is *n individual complex triple product from one 
record in a given lead and BC(fi,f2) is the average complex 
triple product for that same lead: 

bCi(fi,f2) - Xi(fi) * Yi(f2) * ri(fl+f2> 
where l r i(fi+f2) is the complex conjugate of Yi(fi+f2), and 

128 

BC<fi,f 2 ) - jjg £bCi<fi,f 2 ) 

i = 1 

The average real triple product is computed in step 
814 by using the following equations where bri(fi,f2) is an 
individual real triple product from one record in a given 
lead and BR(fi,f2) is the average real triple product for 
that same lead: 

bri(fi,f 2 ) » Pxi(fl) * Pyi(f2) * Pyi(fl+f2) 

128 

BR(fi,f 2 ) » X bri < fl ' f2 > 

i m 1 

In step 816/ the array of auto/crossbispectral density 
values (BD(fi / f2)) is computed using the following equation: 
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BD(fi,f 2 ) « |BC(£i,£ 2 > I 

In step 818, the system computes the array of 
auto/crossbiphase values «p(fi / £ 2 )) using the following 
equation : 

<p(*i/*2) - tan" 1 [Im(BC(fi / f 2 ))/Re(BC(fi,f-))] 
0 < 9 < 2jc (radians) 

In step 820, the system computes the array of 

auto/crossbicoherence values <R(fi,f 2 )) using the following 
equation : 

R(fl,f2) = BO(fi,f2)/[BR<fi,f 2 )] 1/2 
0 < R < 1 

In step 822, the system returns the requested 

auto/cross bispectral density, bicoherence, biphase arrays 

to the monitor module 402. 

Now turning to Fig. 9, a parametric based method for 
producing the autobi spectrum and the crossbispectrum will 
now be described. In steps 902, 904, and 906 the system 
sets the symmetries and time series records in the same man- 
ner as described above in steps 802, 804, and 806 respec- 
tively. The power spectra of Xi(t) and Yi(t) are estimated 
in steps 908, 910, and 912. This estimation method includes 
two major stages, the Autoregressive (AR) model order selec- 
tion and the power spectrum computation for Xi(t) and 
Yi(t). In step 908, the system computes two sequences of 
autocorrelations, {R 2 x(m)} and {R2Y(m)> using the following 
equation . 
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M N-\m[ 

R 2z( m ) = X £ zi(t)zi(t+m), z = X, Y, and m = 0, 1, ... 

i=l t=0 

where M is the number of records of each lead (128 in our 
case) , and N is the number of samples per record (512 in 
our case) , and L is much greater than the possible AR fil- 
ter order (we choose 50) . 

The Final Prediction Errors, FPZx(m) and FPEy(m) are 

calculated for all orders, m » 1, 2, L, by performing a 

Levinson recursion function on each autocorrelation se- 
quence in step 910 in order to find the order of the AR fil 
ter. The locations, Qx and Qy, of the minimum of FPEx(m) 
and FPEy(m) respectively are chosen to be the orders of the 
AR filters of power spectra of Xi(t) and Yi(t) respective- 
ly, i.e., 

FPEx(Qx) = min {fpeM} and FPEy(Qy) = min {FPEy{m)) 

Once the orders of the AR filters for power spectra 
are chosen, the autocorrelation sequences, {R2x(m) } and 
{R2Y(m)}, are entered into Levinson recursion with order Qx 
and Qy, respectively, instead of L. The coefficients, {Cix, 
1=0, 1, . .,,Qx} and {ciy, i - 0,1, ... ,Qy), obtained from 
the recursion are the coefficients of the AR filters for 
power spectra of Xi (t) and Yi(t) respectively. Then, in 
step 912, the power spectra Px(f) and Py(f) are computed as 

the prediction error (o|) divided by square of the mag- 
nitude of the Fourier transform of the coefficients, i.e., 
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Pz(f) = 



Qz 

1 + I c i2 e 
i = 1 



z = X, Y. 



r&*fi .2 



The system estimates the auto/cross bispectrum in 
steps 914, 916, and 918. The estimation process includes 
two major stages: the order selection and bispectrum com- 
putation. In step 914, two sequences of third-order mo- 
ments, {R3x(x)} and {R3y(t)> are computed using the 
following equation. 



M 



52 



* 3z(t) = M~^N 2 X zi(t)z!(t+T), z = X, Y, and x = -L, L 



i=l t-si 



where si » max (1,1-x), s 2 - min (N, N-x) , and L is much 
greater than the possible AR filter orders (e.g. 50). 



In step 916, two super matrices Tx and Ty are formed 
as follows. 



T x = 



K R3zi-L) R3z(-IA-1) „. =? 3 *(0) ^ 
i?3 2 (-L-D R3z(-L) ... Rzzi-l) 

R3z(-2L) R 3 z(-2L+1) .„ =? 3 z<-L) 



z 



X, Y. 



Prom the assumption we made about the AR filter of 
bispectrum, the orders Ox and Oy of the AR filters of 
bispectra of Xi(t) and Yj.(t) are the ranks of the super 
matrices Tx and Ty. Therefore, Ox and Oy are chosen by 
using singular value decomposition. Having found the or- 
ders, we obtain the coefficients of the AR filters of 
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bi spectra by solving the following linear system of equa- 
tions : 



Wo) 


*3z(l) 


... R3z(0z) > 


(l \ 


(*') 




H 32 (0) 


••• R3z(0z-1) 


biz 


0 

r 




R3z(-0z+l) 


... R3z(0) 


bo z z 


0 


V 




J 


v J 


V ) 



where the skewness (J3 Z ) and the coefficients (bi z , 
bo z z) , z ■ X, Y, can be obtained by solving the linear sys- 
tem of equations. 

The auto/cross bispectrum of Xi(t) and Yi(t) are com- 
puted in step 918 as the cubic root of the triple product 

of the skewnesses (Pxpi£y) V3 divided by the triple product of 
the Fourier transforms of the AR filter coefficients 
<H 2 (f)), i.e., 

BC(fi,f 2 ) » ($x&$y) V V HxHf l)HY(f2)Hkf i+f 2 ) 

Hz if ) = 1 + X b± 2 e- j2lcfl , . z = X, Y. 
1 = 1 



and BR(fi,f2) is the real triple product for that same lead: 

BR(fi,f2) = Px(fl) * ?Y(f2) * PY(fl+f2) 

After obtaining power spectrum and auto/cross 
bispectrum, the system computes the bispectral density 
array , the biphase, and the bicoherence in step 920 the 
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same way as in steps 816, 818, 820. In step 922, the system 
returns to the monitor module 402 the requested bispectral 
density, biphase, and bicoherence arrays. 

For illustration purposes Fig. 10 contains sample 
autobispectral arrays showing frequency pairs 0 < fi < 128 
Hz, and 0 < f2 < 64 Hz. A bispectral density array is shown 
in Fig. 10(a) where the Z axis represents the magnitude in 
decibels (db) of the coupled interaction between all ap- 
propriate frequency pairs fi and f2. Recall that the fre- 
quency pairing scheme must adhere to the symmetry rule: 

fi + f2 < N/2 

where N = 256 Hz in this case. A bicoherence array is 
shown in Fig. 10(c) where the Z axis represents the normal- 
ized magnitude in percent (%) of the coupled interaction be- 
tween all appropriate frequency pairs fi and f2. A biphase 
array is shown in Fig. 10(b) where the Z axis represents 
the phase in radians of the coupled interaction between all 
appropriate frequency pairs fi and f2. 

Referring now to Fig 11. a more detailed description 
of the diagnostic index generation module 410 will now be 
provided. In step 1102, the system identifies the type of 
diagnostic assessment in progress. In a preferred embodi- 
ment the four possible options are: 

1. Depth of anesthesia/ pain & surgical stress. 

2. Cerebral ischemia. 

3. Cerebral intoxication (alcohol, narcotics). 

4. Cognitive process evaluation. 
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In step 1104, the system identifies the type of 
bispectral array to pass for use in the diagnostic index 
computation after a user selects a specific lead and array 
type as described above with respect to the user interface 
module 404. There are three (3) possible options for each 
unipolar lead: autobi spectral density; autobiphase; 
autobicoherence. There are also three (3) possible options 
for each set of bipolar leads: crossbi spectral density; 
crossbiphase; crossbicoherence . Since there are 57 (3 X 19 
leads) different types of autobispectral and 9 (3 X 3 sets) 
types of crossbispectral arrays for each one of the 4 diag- 
nostics, the total number of databases is 264. 

In step 1106, the appropriate reference array is 
retrieved from resident memory (or from disk) . Each refer- 
ence array will contain the locations of the frequency 
pairs which are most sensitive to the assessment in 
progress (the generation of the reference arrays and the 
selection of defaults will be discussed later) . In step 
1108, the system adds all data points in the bispectral 
array at the locations identified by the retrieved refer- 
ence array. A counter (NP) of the total number of points 
added is kept. In step 1110, the sum of the data points is 
divided by NP to obtain the single number diagnostic index. 
In step 1112, the program returns to the monitor module 402. 

The predetermined clinical reference arrays referred 
to above are critical to the device's ability to achieve 
clinically relevant diagnostic efficacy. In the following 
section we discuss the process adopted for generating these 
clinical reference arrays. Since a total number of 276 pos- 
sible reference arrays exist, only one will be discussed in 
detail. All other reference arrays are acquired in a 
similar fashion. For illustration purposes the generation 
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of the autobispectral density reference array for monitor- 
ing depth of anesthesia with lead T3 will be reviewed. 

In a first study EEG potentials from a small group of 
medically healthy surgical patients (N) with no known 
neurological disorders are recorded during routine surgery. 
The acquisition procedure described previously 'is followed, 
with the following exception: 

-Band pass filter 0.1 - 500 Hz 

For all patients, two minutes of artifact free EEG 
data are acquired under each of the following conditions: 

- Pre-operative; awake ("control") 

- Deep anesthesia; defined by conventional clini- 
cal standards (intervention or disease state) 

- Post-operative; alert in the recovery room 
(recovery from intervention, or after treatment of disease 
state) 

An autobi spectral density array is generated for lead 
T3 from each one of the three recordings for all patients, 
yielding a total of 3N arrays. The arrays are grouped in 3 
sets of H arrays. The first representing the control state, 
the second representing intervention, and the third repre- 
senting recovery. 

A paired Student's t test is performed on each of 
16,512 data points, comparing the first and second array. 
The resulting 16,512 t values are stored in a two dimen- 
sional array identical in structure to that of the 
bispectral density array. A second paired Student's t test 
is carried out on each of the 16,512 data points, comparing 
the second and third arrays. The resulting 16,512 t values 
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are stored in a second two dimensional array identical in 
structure to that of the bispectral density array. 

All t values not meeting a specific significance test 
or a specific confidence interval in either array are set 
to 0 . In the preferred embodiment all locations with a t 
value not corresponding to a p < 0.0001 are set to 0. Each 
t value from the first t array (Tl(fi,f2)) is compared with 
its corresponding t (T2(fi,f2)) from the second t array. 
One of the following conditions must be met: 
Tl(fi,f 2 ) < 0 < T2(fi,f 2 ) 

or T2(fi,f2) < 0 < Tl(fi / f2) 

If neither one the two conditions is met at a particular 
frequency pair fi,f2 then Tl(fi,f2) * 0 and T2(fi,f2) « 0. 

The application of the above conditions has the effect 
of identifying all of the frequency pair locations that 
change significantly by showing a consistent increase in 
bispectral density value with anesthesia followed by a 
decrease witlf recovery, or a decrease with anesthesia fol- 
lowed by an increase with recovery. 

Finally, the absolute values of the t values in each 
fi/f2 location from the first t array are added to their 
counterpart in the second t array to form a third t array. 
The third t array is an average of the first two and can be 
visually inspected for highly sensitive regions. 

The last step involves sorting the third t array for 
the most sensitive ensemble of frequency pair locations . In 
the preferred embodiment this would consist of the top 25% 
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of all t values. The locations fi,f2 of the most sig- 
nificant t values meeting all of the Above conditions are 
stored in resident memory (or on disk) as one of the 
predetermined reference arrays. This reference array will 
be accessed by the diagnostic index derivation module 410, 
for autobispectral density diagnostic index generation 
during anesthesia/surgery for the location probed by lead 
T3. 

For any particular diagnostic task and any particular 
lead there are 6 possible bispectral arrays 
(autobicoherence, autobispectral density, autobiphase, 
crossbicoherence, cross bispectral density, and 
crossbiphase) which could be examined for diagnostic poten- 
cy. To rank order the reference arrays with respect to 
diagnostic efficacy a second prospective study is con- 
ducted. The conditions under which the study is conducted 
are identical to those of the first except that: a) the fre- 
quency pair locations of interest have already been iden- 
tified and are now followed prospectively and b) the size 
of the study group is now sufficiently large so that sample 
variation of bispectral arrays more closely approximates 
the true variance within the population undergoing the in- 
tervention or suffering from the disease. 

Thus for the example of anesthesia monitoring the EEG 
recording starts during the awake/control state and con- 
tinues uninterrupted through the end of recovery. Con- 
tinuous surgical notes are maintained throughout the 
operation . 

After the completion of the study, continuous diagnos- 
tic indices are generated for the leads of interest for 
each of auto or cross bispectral density, biphase and 
bicoherence arrays. The continuous trends are annotated 
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with the intraoperative notes. A sufficiently large group 
of prospective patients (determined by a statistical power 
test) is used to determine which continuous diagnostic 
index exhibits the greatest diagnostic efficacy on clinical 
grounds. The particular bispectral array used to generate 
this best diagnostic index during a particular diagnostic 
procedure is programmed into the system as the default 
array for diagnostic assessment. 

The following non-limiting example is provided solely 
for illustrative purposes. Twenty (20) patients undergoing 
elective surgery for a variety of orthopedic and 
gynecologic conditions were studied. Standard EEG leads 
were placed in 16 locations according to the International 
10/20 system. Raw EEG signals were acquired, band-pass fil- 
tered (0.1 - 110 Hz) and digitized at a sampling rate of 
256 Hz. EEG recordings were obtained from all patients 
prior to the induction of anesthesia. Patients were then 
anesthetized using standard techniques with a variety of 
anesthetic agents. Continuous EEG recordings were obtained 
during the period of anesthesia induction until the patient 
was judged to be adequately anesthetized for surgery by 
clinical assessment. Intermittent EEG recordings were then 
obtained during the course of the operation. During the 
period of recovery from anesthesia another continuous EEG 
recording was taken. A final recording was obtained when 
the patient was deemed to be "awake" in the recovery room. 
Detailed clinical intra-operative notes of patient status 
were maintained during all phases of EEG recording for sub- 
sequent correlation with bispectral parameters. 

In 10 patients the entire available frequency spectrum 
(0.1 to 110 HZ) was examined for statistically significant 
changes in autobispectral density values from the awake 
state to the deeply anesthetized state and back to the 
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awake state. Pigs. 12 (a) -12(c) show average bispectral den- 
sity arrays (from 10 patients) for each of the three states 
of consciousness. The method for determining statistical 
significance was as outlined above. Pigs. 13 (a) -13(b) show 
the statistical arrays generated by the technique of the 
present invention: the average t array for these 10 sub- 
jects for lead T3 for locations corresponding to a p < 0.05 
(t > 2.26) (Pig. 13(a)) and the average t array for loca- 
tions corresponding to a p < 0.0000003 (t > 10.0) (Fig. 
13 (b) ) . Each t value in the array reflects the consistency 
of change in a bispectral density value through the three 
states for one frequency pair location across all 10 
patients. It is worth noting that virtually no data points 
are significant with a p < 0.0000003 (t > 10.0) in the fre- 
quency pair band of Fl below about 24 Hz and F2 below about 
2 Hz in Pig. 13(b). On the other hand 7,168 locations were 
found to change with a p < 0.0000003 (t > 10.0) in the fre- 
quency pair band of Fl above about 24 Hz and P2 above about 
2 Hz. The top 25% (7,168/4 - 1,792 points) most sig- 
nificant high frequency locations were used to define the 
reference array. The autobi spectral density index was cal- 
culated for each subject from the points defined by the ref- 
erence array as described above in the detailed description 
of the invention. This autobispectral density index was 
then calculated as a continuous function for the continuous 
EEG recordings to assess its behavior during induction, 
intra-operatively, and during recovery. The correlation 
with clinical events during the operation was noted. 

A sample annotated continuous autobispectral density 
diagnostic index for lead T3 during surgery in a prospec- 
tive subject is shown in Fig. 14. The index varies between 
30 and 5 decibels and is quite sensitive to the patient's 
state of consciousness and the onset of painful stimuli. 
Specifically, the index drops with the induction of anes- 
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thetic agents (pentothal and ethrane) to the patient, and 
the index level rises as the patient's leg is being propped 
for surgery. In addition, the index approaches its highest 
value when the patient is awake in the recovery room and 
most likely experiencing post-operative stress. (The gaps 
seen in the index plot correspond to time periods when EEG 
recordings were not being taken) . 

Similarly the above analytic process is used to 
generate the reference databases for cerebral ischemia, 
degrees of intoxication and normal or abnormal cerebral 
processes. In quantitatively detecting any of these 
cerebral phenomina, the system compares a number of 
autobispectral and crossbispectral EEG data from subjects 
in the normal state to clinically identified extremes of a 
certain physiologic state (awake vs anesthetized, sober vs 
intoxicated, perfused vs ischemic, at rest mentally vs 
thinking, normal vs retarded, etc..)* The comparison util- 
izes a statistical approach to identify the bispectral data 
points that are^most sensitive to the particular 
physiologic state in question. The frequency pair loca- 
tions of the most sensitive data points are identified and 
stored in a database for reference purposes. When a diag- 
nosis is to be carried out, the average of all the data 
points defined by the reference array is obtained for the 
subject undergoing the study. This average is used as a 
diagnostic index and is compared to a list of indices char- 
acteristic of each state by the operator or the system. 

In addition to quantifying the depth and adequacy of 
anesthesia, pain responses during surgical stress, acute 
and chronic cerebral ischemia, level of consciousness, de- 
gree of cerebral intoxication and normal or abnormal cogni- 
tive processes, the system and method of the present 
invention may also be used to assess a myriad of cerebral 
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phenomena based on the acquisition and processing of EEG 
signals into various bispactral arrays which are then com- 
pared to appropriate reference arrays. 

Although bi spectral analytic techniques in the frequen- 
cy domain have been applied to the EEG signal, as was dis- 
cussed in the Background above, parametric approaches to 
the estimation of bispectral values have not. Furthermore 
no bispectral technique has ever been demonstrated to be 
useful for any diagnostic purpose. Other techniques for the 
quantification of the depth of anesthesia or the detection 
and quantification of cerebral ischemia intraoperatively 
remain qualitative and limited in their overall utility and 
acceptance in practice. Specifically, the system and method 
of the present invention examines various bispectral values 
across all frequency pairs in a frequency range hitherto ig- 
nored by those knowledgeable in the art and uses the summed 
degree of changes in autobicoherence/autobispectral den- 
sity/autobiphase, crossbicoherence/crossbispectral den- 
sity/crossbiphase at a limited number of frequency 
locations as an index of physiological perturbation. The 
system and method utilize various bispectral arrays of 
defined clinical populations to define the locations of the 
subset of frequencies used to calculate this index. Refer- 
ence clinical arrays are further utilized to assess the 
meaning of this index and to measure the significance of 
deviations of this index from normality. This allows the 
quantitative gauging of the disturbances in cerebral func- 
tion, whether due to anesthesia, intoxicants or ischemia 
for any particular EEG lead position. The system and method 
disclosed herein also define the graphic display of the 
diagnostic index, whether on graphics screen or on paper, 
whether in real-time or in digital archive. 
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While the foregoing invention has been described with 
reference to its preferred embodiments/ various alterations 
and modifications will occur to those skilled in the art. 
All such alterations and modifications are intended to fall 
within the scope of the appended claims. 
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What is claimed is: 

1. A method of noninvasively detecting cerebral 
phenomena comprising the steps of: 

acquiring elect roencephalographic signals through 
at least one electrode from a body surface of a subject 
being analyzed; 

filtering said electroencephalograph^ signals to 
obtain filtered signals having frequencies between 2 and 
500 hertz; 

dividing said filtered signals into a plurality 
of equally sized data records; 

characterizing dynamic phase relations within 
said filtered signals by processing said filtered signals 
to generate bispectral values; 

comparing said generated bispectral values to ref- 
erence values to derive a diagnostic index that quantifies 
the detected cerebral phenomena. 

2. The method of noninvasively detecting cerebral 
phenomena of claim 1 whferein the step of acquiring 
electroencephalographic signals further comprises the step 
of attaching electrodes to a head of the subject being 
analyzed in order to obtain a unipolar electroen- 
cephalographic signal from each region of interest of both 
left and right hemispheres of the subject's brain to which 
said electrodes are attached. 

3. The method of noninvasively detecting cerebral 
phenomena of claim 1 wherein said bispectral values 
generated in said step of characterizing said dynamic phase 
relations are autobispectral density values. 



WO 90/11718 



PCT/US90/01378 



-38- 

4. The method of noninvasively detecting cerebral 
phenomena of claim 1 wherein said bi spectral values 
generated in said step of characterizing said dynamic phase 
relations are autobicoherence values. 

5* The method of noninvasively detecting cerebral 
phenomena of claim 1 wherein said bispectral values 
generated in said step of characterizing said dynamic phase 
relations are autobiphase values. 

6. The method of noninvasively detecting cerebral 
phenomena of claim 1 wherein said step of acquiring 
electroencephalograph^ signals further comprises the step 
of attaching electrodes to the head of the subject being 
analyzed in order to obtain bipolar data sets of electroen- 
cephalograph^ signals from left and right hemispheres of 
the subject's brain to which said electrodes are attached. 

7. The method of noninvasively detecting cerebral 
phenomena of claim 6 wherein one bipolar data set is ac- 
quired from a frontal left hemisphere of the subject's 
brain and another bipolar data set is acquired from a fron- 
tal right hemisphere of the subject's brain. 

8. The method of noninvasively detecting cerebral 
phenomena of claim 6 wherein one bipolar data set is ac- 
quired from a left occipital region of the subject's brain 
and another bipolar data set is acquired from a right oc- 
cipital regioi of the subject's brain. 
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9. The method of noninvasively detecting cerebral 
phenomena of claim 6 wherein one bipolar data set is ac- 
quired from a left parietal region of the subject's brain 
and another bipolar data set is acquired from a right par 
tal region of the subject's brain. 



10. The method of noninvasively detecting cerebral 
phenomena of claim 1 wherein said bispectral values 
generated in said step of characterising said dynamic phase 
relations are crossbispectral density values. 

11. The method of noninvasively detecting 
cerebral phenomena of claim 1 wherein said bispectral 
values generated in said step of characterizing said 
dynamic phase relations are crossbicoherence values. 

12. The method of noninvasively detecting 
cerebral phenomena of claim 1 wherein said bispectral 
values are generated in said step of characterizing said 
dynamic phase relations are crossbiphase values. 
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13. The method of noninvasively detecting 
cerebral phenomena of claim 3 where said step of generating 
autobispectrtfl -density values comprises the steps of: 

computing fast Fourier transforms Xi(f) and Yi(f) 
of each of said data records i; 

computing power spectra Pxi(f) and Pyi(f) of said 
data records i by squaring the magnitude of each element of 
said fast Fourier transforms Xi(f) and Yi(f) respectively; 

computing for at least one electrode an average 
complex triple product of data records acquired by said at 
least one electrode; 

computing for said at least one electrode an 
average real triple product of data records acquired by 
said at least one electrode; 

computing for said at least one electrode an 
autobi spectral density value as the absolute value of said 
average complex triple product for said electrode. 

14. The method of noninvasively detecting cerebral 
phenomena of claim 13 further comprising the step of comput- 
ing an autobiphase value <p(fi,f2) for at least one 
electrode such that: 

9(*i/*2) - tan" 1 [Im(BC(fi / f2))/Re(BC(fi,f2))] 
where BC(fi / f2) is the average complex triple product for 
an electrode, and f i and f2 designate limits of the frequen- 
cy range over which the autobiphase computation is carried 
out. 
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15. The method of noninvasively detecting cerebral 
phenomena o£ claim 13 further comprising the step of comput- 
ing an autobicoherence value R(fi,f2> for at least one 
electrode such that 

R(fl,f2) = BD(t lf t Z )/[BR(f lf f 2 )] 1/2 
where BD(fi / f2) is the aut obi spectral density value for 
said electrode, BR(fi,f 2 ) is the average real triple 
product for the same electrode, and fi and f2 designate 
limits of the frequency range over which the 
autobicoherence computation is carried out. 

16. The method of noninvasively detecting cerebral 
phenomena of claim 10 where said step of generating said 
crossbispectral density values comprises the steps of: 

computing fast Fourier transforms Xi(f) and Yi(f) 
of said data records i; 

computing power spectra Pxi(f) and PYi(f) of said 
data records by squaring the magnitude of elements of said 
fast Fourier transforms Xi(f) and Yi(f) respectively; 

computing for at least one electrode pair an 
average complex triple product of all data records acquired 
by said at least one electrode pair; 

computing for said at least one electrode pair an 
average real triple product of all data records acquired by 
each of said at least one electrode pair; 

computing for said at least one electrode pair a 
crossbispectral density value as the absolute value of the 
average complex triple product for said electrode pair. 
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17. The method of noninvasively detecting cerebral 
phenomena of claim 16 further comprising the step of comput- 
ing a crossbiphase value (p(fi,f2) for said at least one 
electrode pair such that: 

<p(fl,f2> * tan" 1 [Im(BC(fi / f 2 ))/Re(BC(fi / f 2 ))] 
where BC(fi,f2> is the average complex triple product for 
an electrode pair, and fi and f2 designate limits of the 
frequency range over which the crossbiphase computation is 
carried out. 

18. The method of noninvasively detecting cerebral 
phenomena of claim 17 further comprising the step of comput- 
ing a crossbicoherence value R<fi,f2) for said at least one 
electrode pair such that 

R(fl,f2) - BD<fi,f 2 ) /[BR(fi,f 2 )] 1/2 
where BD(fi,f2) is the ctossbi spectral density value for an 
electrode pair, BR(fi / f2) is the average real triple 
product for the same electrode pair, and f i and f 2 desig- 
nate limits of the frequency range over which 
crossbicoherence computation is carried out. 
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19. The method of noninvasively detecting cerebral 
phenomena of claim 3 wherein said step of generating 
aut obi spectral density values comprises the steps of: 

computing autocorrelation sequences R2x(m) and 
R2Y(m) of all data records acquired by at least one 
electrode; 

determining the orders and coefficients of~ 
parametric models for power spectra of data records ac- 
quired by said at least one electrode; 

computing power spectra Px(f) and Pyff) of data 
records acquired by said at least one electrode; 

computing third order moment sequences R3x(x) and 
R3y(t) of data records acquired by said at least one 
electrode; 

determining the orders and coefficients of 
parametric models of the bispectra of data records acquired 
by said at least one electrode; 

computing for said at least one electrode a 
bispectrum of data records acquired by said at least one 
electrode . 

i 

20. The method of noninvasively detecting cerebral 
phenomena of claim 19 wherein said bispectrum is 
autobi spectrum and further comprising the step of computing 
an autobispectral density value for at least one electrode 
as the absolute value of the bispectrum of all data records 
for said electrode. 
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21. The method of noninvasively detecting cerebral 
phenomena of claim 19 further comprising the step of comput- 
ing an autobiphase value <p(fi,f2) for said at least one 
electrode such that: 

<p(fl,f2) ■ tan" 1 [Im(BC(fi / f 2 ))/Re(BC(fi / f2))] 
where BC(fi,f2) is the bispectrum for an electrode, and fi 
and ±2 designate limits of the frequency range over which 
the bispectral computation is carried out. 

22. The method of noninvasively detecting cerebral 
phenomena of claim 20 further comprising the steps of: 

computing for at least one electrode a real 
triple product of data records acquired by said at least 
one electrode; 

computing an autobi coherence value R(fi,f2) for 
said at least one electrode such that 

R(fl,f2> « BD(fi / f2)/[BR(fi / f 2 )] 1/2 

where BD(fi,f2) is the autobi spectral density value for an 
electrode, BRffi,f 2 ) is the real triple product for the 
same electrode, and f i and f 2 designate limits of the fre- 
quency range over which bispectral computation is carried 
out. 

23. The method of noninvasively detecting cerebral 
phenomena of claim 19 wherein said bispectrum is 
crossbi spectrum and further comprising the step of comput- 
ing a crossbi spectral density value for each electrode pair 
as the absolute value of the crossbispectrum of all data 
records for each said electrode pair. 
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24. The method of noninvasively detecting cerebral 
phenomena of claim 23 further comprising the step of comput- 
ing a crossbiphase value <p<*i/f2) for each of said at 
least one electrode pair such that: 

<P(fi,f2) = tan" 1 [Im(BC(fi,f 2 ))/Re(BC(fi,f2))] 
where BCCfi,f2) is the crossbispectrum for an pair, and fi 
and f2 designate limits of the frequency range over which 
the bispectral computation is carried out. 

25. The method of noninvasively detecting cerebral 
phenomena of claim 24 further comprising the step of comput- 
ing a crossbicoherence value R(fi,f2) for each of said at 
least one electrode pair such that 

R(fl,f2) « BD(fi,f2)/[BR(f 1/ f 2 )] 1/2 
where BD(fi,f 2 ) is the crossbispectral density value for an 
electrode pair, BR(fi,f 2 ) is the real triple product for 
the same electrode pair, and fi and f2 designate limits of 
the frequency range over which bispectral computation is 
carried out. 

26. The method for noninvasively detecting cerebral 
phenomena of claim 1 wherein said step of acquiring 
electroencephalographic signals further comprises the step 
of analyzing said signals to determine lead failure. 
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27 • ^The method of noninvasively detecting cerebral 
phenomena of claim 1 wherein said step of comparing further 
comprises the steps of: 

organizing said generated bispectral values in at 
least one array of bispectral values ; 

selecting a physical phenomena to be diagnosed; 

retrieving an appropriate bispectral reference 
array from a resident memory, said reference array contain- 
ing frequency pairs that are most sensitive to the physical 
phenomena to be diagnosed; 

adding data values in locations of each of said 
at least one array of bispectral values that are identified 
by the retrieved reference array as being locations contain- 
ing data of significance to obtain a sum of said sig- 
nificant locations; 

averaging the values stored in said significant 
locations to generate a diagnostic index relating to the 
cerebral phenomena to be detected. 
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28. The method of noniavasively detecting cerebral 
phenomeaa of claim 1 further comprising the steps of: 

geaeratiag three arrays of bispectral data for 
each of three different states of the subject; 

performing a paired Student's t test comparing 
data in a first and a second array of said three arrays of 
bispectral data to produce a first t array and performing a 
paired Student's t test comparing data in said second and a 
third array of said three arrays of bispectral data to 
produce a second t array; 

comparing data values in said first t array with 
data values in corresponding locations in said second t 
array; 

identifying those corresponding locations in said 
first and second t arrays that differ by more than a pre- 
selected amount, said identified locations representing 
those locations that are significant for detecting the 
cerebral phenomena. 

29. The method of noninvasively detect iag cerebral 
phenomena of claim 1 wherein the cerebral phenomenon being 
detected is the depth of anesthesia in the subject being 
analyzed. 

30. The method of noninvasively detecting cerebral 
phenomena of claim 1 wherein the cerebral phenomena being 
detected are pain responses during surgical stress in the 
subject being analyzed. 

31. The method of noninvasively detecting cerebral 
phenomena of claim 1 wherein the cerebral phenomenon being 
detected is acute and chronic ischemia or infarction in the 
subject being analyzed. 
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32. The method of noninvasively detecting cerebral 
phenomena of claim 1 wherein the cerebral phenomenon being 
detected is the level of consciousness of the subject being 
analyzed. 

33. The method of noninvasively detecting cerebral 
phenomena of claim 1 wherein the physical phenomenon being 
detected is the degree of cerebral intoxication of the sub- 
ject being analyzed. 

34. The method of noninvasively detecting physical 
phenomena of claim 1 wherein the physical phenomena being 
detected are normal or abnormal cognitive processes. 

35. A system for noninvasively detecting cerebral 
phenomena comprising: 

means for acquiring electroencephalograph^ sig- 
nals through at least one electrode from a body surface of 
a subject being analyzed; 

means for filtering said electroencephalograph^ 
signals to eliminate those signals having frequencies less 
than 2 hertz or frequencies greater than 500 hertz; 

means for dividing said filtered signals into a 
plurality of equally sized data records; 

means for generating bispectral values capable of 
characterizing dynamic phase relations within said filtered 
electroencephalograph^ signals; 

means for comparing said generated bispectral 
values to reference values in order to derive a diagnostic 
index that quantifies the detected cerebral phenomena. 
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36. The system for noninvaaively detecting cerebral 
phenomena of claim 35 further comprising a plurality of 
said means for acquiring electroencephalographic signals, 
each of said means for acquiring electroencephalographic 
signals being connected to said means for filtering. 

37. The system for. noninvaaively detecting cerebral 
phenomena of claim 36 wherein said plurality of said means 
for acquiring electroencephalographic signals is a 
plurality of electrodes attachable to a head of a subject 
being analyzed to obtain a unipolar electroencephalographic 
signal from each of a plurality of regions of interest on 
both left and right hemispheres of the subject's brain. 

38. The system for noninvasively detecting 
cerebral phenomena of claim 35 wherein said means for ac- 
quiring encephalography signals comprises: 

a plurality of surface electrodes for mounting on 
a surface of a head of the subject being analyzed; 

means for providing elactrosurgery protection in- 
cluding a radio frequency filter for limiting radio frequen- 
cy current through said electrodes ; 

means for providing defibrillator protection for 
limiting voltage to said amplifiar during a discharge; 

means for amplifying said filtered signals for a 
high gain in order to maximize the dynamic range for high 
frequency, low energy wave components of said filtered sig- 
nals; 

means for arranging said filtered signals in a 
serial stream and for feeding said signals to an analog-to- 
digital converter to convert said signals to digital sig- 
nals. 
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39. The system for noninvasively detecting cerebral 
phenomena of claim 35 further comprising means for analyz- 
ing said signals received by each electrode in order to 
detect electrode failure. 

40. The system for noninvasively detecting cerebral 
phenomena of claim 35 wherein said bispectral values are 
autobi spectral density values and further comprising means 
for organizing said autobi spectral density values in at 
least one array of autobi spectral density values. 

41. The system for noninvasively detecting cerebral 
phenomena of claim 35 wherein said bispectral values are 
autobicoherence variables and further comprising means for 
organizing said autobicoherence values in at least one 
array of autobicoherence values. 

42. The system for noninvasively detecting cerebral 
phenomena of claim 35 wherein said bispectral values are 
autobiphase values and further comprising means for organiz- 
ing said autobiphase values in at least one array of 
autobiphase values. 

43. The system for noninvasively detecting cerebral 
phenomena of claim 35 wherein said bispectral values are 
crossbispectral density values and further comprising means 
for organizing said crossbispectral density values in at 
least one array of crossbispectral density values. 

44. The system for noninvasively detecting cerebral 
phenomena of claim 35 wherein said bispectral values are 
crossbicoherence values and further comprising means for 
organizing said crossbicoherence values in at least one 
array of autobicoherence values. 
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45. The system for noninvasively detecting cerebral 
phenomena of claim 35 wherein said bispectral variables a 
crossbiphase values and further comprising means for or- 
ganizing said crossbiphase values in at least one array o 
autobiphase values. 



46. The system for noninvasively detecting cerebral 
phenomena of claim 35 wherein said means for acquiring en- 
cephalography signals further comprises means for obtain- 
ing bipolar data sets of electroencephalograph^ signals 
from different regions of a brain of said subject. 

47. The system for noninvasively detecting cerebral 
phenomena of claim 40 where said means for generating at 
least one array of autobispectral density values comprises: 

means for computing fast Pourier transforms Xi(f) 
and Yi(f) of each of said data records i; 

means for computing power spectra Pxi(f) and 
PYi(f) of said data records i by squaring the magnitude of 
elements of said fast Fourier transforms Xi(f) and Yi(f) 
respectively; 

means for computing for said at least one 
electrode an average complex triple product of data records 
acquired by said at least one electrode; 

means for computing for said at least one 
electrode an average real triple product of data records ac- 
quired by said at least one electrode; 

means for computing for said at least one 
electrode an autobispectral density value as the absolute 
value of the average complex triple product for said 
electrode . 
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48. The system of noninvasively detecting cerebral 
phenomena of claim 47 further comprising means for comput- 
ing an aut&biphase value <p(fi,f2) for said at least one 
electrode suci that: 

<p(fl,f2) = tan" 1 [lm(BC (fi,f2))/Re (BC(£i,f 2 ))] 

where BC (f i, f2) is the average complex triple product for 
an electrode, and fi and f2 designate limits of the fre- 
quency range over which the autobiphase computation is car- 
ried out. 

49. The system for noninvasively detecting cerebral 
phenomena of claim 47 further comprising means for comput- 
ing an *ut obi coherence value R(fi,f2) for said at least one 
electrode such that 

R(fl,f2) » BD(£i,f2)/[BR(£i r £ 2 )] 1/2 

where BD (f i, f 2) is the aut obi spectral density value for 
said electrode, BR(fi,f2) is the average real triple 
product for the same electrode, and f 1 and f 2 designate 
limits of the frequency range over which the 
autobicoherence computation is carried out. 
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50. The system for noninvasively detecting cerebral 
phenomena of claim 40 where said means for generating at 
least one array of crossbi spectral density values comprises 
means for computing fast Fourier transforms Xi(f) 
and Yi(f) of each of said data records i; 

means for computing power spectra Pxi(f) and 
PYi(f) of said data records by squaring the magnitude of 
elements of said fast Pourier transforms Xi(f) and Yi(f) 
respectively; 

means for computing for at least one electrode 
pair an average complex triple product of all data records 
acquired by each of said at least one electrode pair; 

means for computing for said at least one 
electrode an average real triple product of all data 
records acquired by for each of said at least one electrode 
pair; 

means for computing for said at least one 
electrode pair a crossbispectral density value as the ab- 
solute value of the average complex triple product for said 
electrode pair. 

< 

51. The system for noninvasively detecting cerebral 
phenomena of claim 50 further comprising means for comput- 
ing a crossbiphase value <p(fi,f 2 ) for said at least one 
electrode pair such that: 

9<*i/ f2) » tan" 1 [Im(BC(f 1 ,f 2 ))/Re(BC(fi / f2))l 
where BC(fi,f 2 ) is the average complex triple product for 
an electrode pair, and fi and f 2 designate limits of the 
frequency range over which the crossbiphase computation is 
carried out. 
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52 . The system for noninvasively detecting cerebral 
phenomena of claim 51 further comprising means for comput- 
ing an crossbicoherence value R(fi/f2) for said at least 
one electrode pair such that 

R(fl,f2> - BD(fi / f2)/[8R(fi / f 2 )] 1/2 

where BD (f i, f2) is the crossbispectral density value for an 
electrode pair, BR(fi,f2) is the average real triple 
product for the same electrode pair, and fi and f2 desig- 
nate limits of the frequency range over which the 
crossbicoherence computation is carried out. 

53. The system for noninvasively detecting cerebral 
phenomena of claim 40 wherein said means for generating at 
least one array of autobispectral density values comprises: 

means for computing autocorrelation sequences 
R2x(m) and R2y (m) of all data records acquired by at least 
one electrode; 

means for determining the orders and coefficients 
of parametric models for power spectra of data records ac- 
quired by said at least one electrode; 

means for computing power spectra Px(f) and Py(f) 
of all data records acquired by said at least one electrode; 

means for computing third order moment sequences 
R3x(T) and R3y(?) of data records acquired by said at least 
one electrode; 

means for determining the orders and coefficients 
of parametric models of the bispectra of data records ac- 
quired by said at least one electrode; 

means for computing for said at least one 
electrode a bispectrum of data records acquired by said at 
least one electrode; 
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54. The system for noninvasively detecting cerebral 
phenomena of claim 53 wherein said bispectrum is 
autobispectrum and further comprising means for computing 
an autobispectral density value for each electrode as the 
absolute value of the bispectrum of data records for said 
electrode . 

55. The system for noninvasively detecting cerebral 
phenomena of claim 53 further comprising means for comput- 
ing an autobiphase value (p(fi,f2) for at least one 
electrode such that: 

<p(fi,f2) = tan" 1 [Im(BC(fi / f 2 ))/Re(BC(fi f f 2 ))] 
where BC(fi,f2) is the bispectrum for an electrode, and fi 
and f2 designate limits of the frequency range over which 
the autobiphase computation is carried out. 

56. The system for noninvasively detecting cerebral 
phenomena of claim 54 further comprising: 

means for computing for at least one electrode a 
real triple product of all data records acquired by said at 
least one electrode; 

means for computing an autobicoherence value 
R(fi,f2) for said at least one electrode such that 

R(fi,f2) = BD<fi,f 2 )/[BR<fi,f 2 )] 1/2 

where BD(fi,f 2 ) is the autobispectral density value for an 
electrode, BR(fi,f2) is the real triple product for the 
same electrode, and fi and f2 designate limits of the fre- 
quency range over which autobicoherence computation is car- 
ried out. 
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57. The system for noninvasively detecting cerebral 
phenomena of claim 53 wherein said bispectrum is 
crossbi spectrum and further comprising means for computing 
a crossbispectral density value for an electrode pair as 
the absolute value of the bispectrum of data records for 
said electrode pair. 

58. The system for noninvasively detecting cerebral 
phenomena of claim 53 further comprising means for comput- 
ing a crossbiphase value <p(£i,£2) for at least one 
electrode pair such that: 

9(fi/f2) = tan" 1 [Im(BC(fi / f 2 ))/Re(BC(fi / f2))3 
where BC(fi,f2) is the crossbi spectrum for an electrode 
pair, and fi and f2 designate limits of the frequency range 
over which the -crossbiphase computation is carried out. 

59. The- system for noninvasively detecting cerebral 
phenomena of claim 57 further comprising: 

means for computing for at least one electrode a 
real triple product of all data records acquired by said at 
least one electrode; 

means for computing a crossbi coherence value 
R(fi/f 2 ) for said at least one electrode pair such that 

R(fl,f2) - BD(fi / f 2 )/[BR(fi / f 2 )] 1/2 
where BD(fi, ti) is the crossbispectral density value for an 
electrode pair/ BR(fi,f2) is the real triple product for 
the same electrode pair, and f i and f2 designate limits of 
the frequency range over which the crossbicoherence computa 
tion is carried out. 
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60. The system for noninvasively detecting cerebral 
phenomena of claim 35 wherein said means for comparing fur- 
ther comprises: 

means for organizing said generated bispectral 
values in an array of bispectral variables; 

means for selecting a- physical phenomena to be 

diagnosed; 

means for retrieving an appropriate bispectral 
reference array from a resident memory, said reference 
array containing frequency pairs that are most sensitive to 
the physical phenomena to be diagnosed; 

means for adding data values in locations of each 
of said at least one array of bispectral values that are 
identified by the retrieved reference array as being loca- 
tions containing data of significance to obtain a sum of 
said significant locations; 

means for averaging the values stored in said sig- 
nificant locations to generate a diagnostic index relating 
to the cerebral phenomena to be detected. 
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61. The system for noninvasively detecting cerebral 
phenomena of claim 1 further comprising: 

means for generating three arrays of bi spectral 
data for each of three different states of the subject; 

means for performing a paired Student' s t test 
comparing the data in a first and a second array of said 
three arrays of bispectral data to -produce a first t array 
and performing a paired Student's t test comparing the data 
in said second and a third array of said three arrays of 
bispectral data to produce a second t array; 

means for comparing each data value in said first 
t array with data values in corresponding locations in said 
second t array ; 

means for identifying those corresponding loca- 
tions in said first and second t arrays that differ by more 
than a pre-selected amount, said identified locations repre- 
senting those locations that are significant for detecting 
the cerebral phenomena. 

62. The system for noninvasively detecting cerebral 
phenomena of claim 35 further comprising means for display- 
ing a representation of a subject's head, divided into a 
selected number of sections, said means for displaying in- 
cluding means for displaying a compressed continuous trac- 
ing of a computer diagnostic index determined from the 
signals acquired from an electrode positioned at a loca- 
tion represented by said section. 

63. The system for noninvasively detecting cerebral 
phenomena of claim 61 wherein each displayed section in- 
cludes a background that is color-coded to reflect possible 
values in a range of a selected diagnostic index. 
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64. The method of noninvasively detecting cerebral 
phenomena of claim 1 further comprising the steps of: 

generating three arrays of bispectral data for 
each of three different states of the subject; 

performing statistical operations on said three 
arrays of bispectral data in order to identify those loca- 
tions in said arrays that are significant for detecting the 
cerebral phenomena. 

65. The system for noninvasively detecting cerebral 
phenomena of claim 1 further comprising: 

means for generating three arrays of bispectral 
data for each of three different states of the subject; 

means for statistically analyzing said arrays of 
bispectral data in order to identify those locations in 
said arrays that are significant for detecting the cerebral 
phenomena . 
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AMENDED CLAIMS 
[received by the International Bureau 
on 19 September 1990 (19.09.90); 
original ^claims 31, 38, 61 and 63 amended; 
new claims 66-70 added; other claims anchanged (5 pages)] 

28. The method of noninvasively detecting 
cerebral phenomena of claim 1 further comprising the 
steps of: 

generating three arrays of bispectral data for 
each of three different states of the subject; 

performing a paired Student's t test comparing 
data in a first and a second array of said three 
arrays of bispectral data to produce a first t array 
and performing a paired student's t test comparing 
data in said second and a third array of said three 
arrays of bispectral data to produce a second t array; 

comparing data values in said first t array with 
data values in corresponding locations in said second 
t array; 

identifying those corresponding locations in said 
first and second t arrays that differ by more than a 
preselected amount, said identified locations 
representing those locations that are significant for 
detecting the cerebral phenomena. 

29. The method of noninvasively detecting 
cerebral phenomena of claim 1 wherein the cerebral 
phenomenon being detected is the depth of anesthesia 
in the subject being analyzed. 

30. The method of noninvasively detecting 
cerebral pheomena of claim 1 wherein the cerebral 
phenomena being detected are pain responses during 
surgical stress in the subject being analyzed; 

31. (Amended) The method of noninvasively 
detecting cerebral phenomena of claim 1 wherein the 
cerebral phenomenon being detected is acute ischemia 
or infarction in the subject being analyzed. 
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36. The system for noninvasively detecting 
cerebral phenomena of claim 35 further comprising a 
plurality of said means for acquiring 
electroencephalograph^ signals, each of said means 
for acquiring electroencephalographic signals being 
connected to said means for filtering. 

37. The system for noninvasively detecting 
cerebral phenomena of claim 36 wherein said plurality 
of said means for acquiring electroencephalographic 
signals is a plurality of electrodes attachable to a 
head of a subject being analyzed to obtain a unipolar 
electroencephalographic signal from each of a 
plurality of regions of interest on both left and 
right hemispheres of the subject's brain. 

38. (Amended) The system for noninvasively 
detecting cerebral phenomena of claim 35 wherein said 
means for acquiring electroencephalographic signals 
comprises : 

a plurality of surface electrodes for mounting on 
a surface of a head of the subject being analyzed; 

means for providing electrosurgery protection 
including a radio frequency filter for limiting radio 
frequency current through said electrodes; 

means for providing defibrillator protection for 
limiting voltage to said amplifier during a discharge; 

means for amplifying said filtered signals for a 
high gain in order to maximize the dynamic range for 
high frequency , low energy wave components of said 
filtered signals; 

means for feeding said signals to an analog-to- 
digitial converter to convert said signals to digital 
signals . 
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61. (Amended) The system for noninvasively 
detecting cerebral pheomena of claim 35 further 
comprising : 

means for generating three arrays of bispectral 
data for each of three different states of the 
subject; 

means for performing a paired Student's t test 
comparing the data in a first and a second array of 
said three arrays of bispectral data to produce a 
first t array and performing a paired Student's t test 
comparing the data in said second and a third array of 
said three arrays of bispectral data to produce a 
second t array; 

means for comparing each data value in said first 
t array with data values in corresponding locations in 
said second t array; 

means for identifying those corresponding 
locations in said first and second t arrays that 
differ by more than a pre-selected amount, said 
identified locations representing those locations that 
are significant for detecting the cerebral pheomena, 

62. The system for noninvasively detecting 
cerebral phenomena of claim 35 further comprising 
means for displaying a representation of a subject's 
head, divided into a selected number of sections, said 
means for displaying including means for displaying a 
compressed continuous tracing of a computer diagnostic 
index determined from the signals acquired from an 
electrode positioned at a location represented by said 
section. 
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63. (Amended) The system for noninvasively 
detecting cerebral phenomena of claim 61 wherein each 
displayed section includes a background of one of a 
plurality of colors, each of which colors is unique to 
a distinct selected range of possible values of a 
selected diagnostic index. 
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66. (new) The method of noninvasively detecting 
cerebral phenomena of claim 1 wherein the cerebral 
phenomenon being detected is chronic ischemia or 
infarction in the subject being analyzed. 

67. (new) The method of noninvasively detecting 
cerebral phenomena of claim 1 wherein said bispectral 
values are generated by computing the Fourier 
transform of the third order autocorreltaiton function 
of said filtered signals. 

68. (new) The method of noninvasively detecting 
cerebral pheonemona of claim 1 wherein said bispectral 
values are generated by computing the Fourier 
transform of the third order crosscorrelation function 
of said filtered signals. 

69. (new) The system for noninvasively 
detecting cerebral phenomena of claim 35 wherein said 
means for generating bispectral values comprises a 
means for computing the Fourier transform of the third 
order autocorrelation function of said filtered 
signals . 

70. (new) The system for noninvasively 
detecting cerebral phenomena of claim 35 wherein said 
means for generating bispectral values comprises a 
means for computing the Fourier transform of the third 
order crosscorrelation function of said filtered 
signals. 
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STATEMENT UNDER ARTICLE 19 



The amendments to claims 31, 38, 61 and 63 were 
made to correct obvious dependence errors in the 
claims and to remove unneeded language, and support 
for the new claims 66 through 70 can be found 
throughout the applications as originally filed as 
bispectral analysis is defined as the Fourier 
transform of the third order autocorrelation function. 
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